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1. Introduction 



We consider the one-dimensional Fokker-Planck equation 
d V . . „ d 



dx'^ dx 
or the equivalent Schrodinger equation 

d^ 



<t>{x) + 2—[f{x)<t>{x)]^e<l>{x) (1.1) 



^^2V'(x) + Vs{x)iP{x) = k'iPix). (1.2) 

Equation (1.1) describes the diffusion of particles in an external potential V{x), from 
which the function f{x) in (1.1) is defined by 

m = (1-3) 

The Schrodinger potential Vs and the function -0 of (1.2) are are related to / and 4> by 
Vs{x)^f\x) + f{x) (1.4) 

and 

0(x) = e-^(^)/V(a:). (1.5) 

We shall always assume that Im /c > 0. We define the Green function Gs{x, y; k) for the 
Schrodinger equation as the function satisfying 



-Vs{x) + k' 



Gs{x,y;k)^5{x-y) (1.6) 



dx'^ 

with the boundary conditions Gs{x,y;k) — >■ as \x — y\ — >■ oo for Im/c > 0. For 
Im A; = 0, we define Gs{x, y; k) = lim^o Gs{x, y;k + ie). Without loss of generality, we 
may suppose that x > y. 

In a series of previous papers [1,2], we discussed a method for calculating the 
expansion of Gs{k) in powers of k. In [1], we studied the cases in which the potential 
V{x) cither converges to a finite limit of diverges to infinity as a; — )■ ±cxd. In [2], we 
dealt with periodic potentials satisfying V{x + L) = V{x). In this paper, we shall deal 
with asymptotically periodic potentials, i.e. potentials V{x) that approach a periodic 
function as a; — )■ ±oo. In solid state physics, impurities in a crystal are described by 
this type of potentials. The study of asymptotically periodic potentials is important for 
the application in physics, and there is a fair amount of literature on this subject [3-12]. 
However, there has not yet been a systematic analysis of the the low-energy behavior of 
the Green function up to high orders in k. In this paper, we shall show that the method 
introduced in [1] and [2] can be extended to the asymptotically periodic case, enabling 
us to obtain the expansion of the Green function up to any order in k. 

We assume that the potential V{x) is a real- valued function which is piecewise 
continuously differentiable. (Note that V{x) may have have jump discontinuities. See 
footnote 1 of [2] and footnote 1 of [1].) We also assume that V can be expressed as a 
sum of two functions: 

V{x) = Vp{x) + Va{x), (1.7) 
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where T/p is a periodic function satisfying 

Vp(x + L)^V^(x), (1.8a) 
and Va is a function such that 

hm VAix) = 0. (1.86) 

x— >±oo 

Corresponding to (1.7), we assume that the function / (equation (1.3)) can be written 
as 

/(x) = /p(x)+/a(x), (1.9) 

where fp{x + L) = fp{x) and \im^^±^ /a(x) = 0. 

As we will see, the Green function can be expanded in terms of k as 

Gs{x, y; k) = (iA;)~^5r_i(a:, y) + go{x, y) + \kgi{x, y) + {\kfg2{x, y) 

+ --- + (ifc)^(/^(a;,|/) + o(A:^) (1.10) 

if V/^ e L]v, where L\ denotes the set of functions v{x) satisfying 

{l + \x\^)\v{x)\dx <oo. (1.11) 



/oo 
-oo 



In this paper, we shall discuss the method for systematically calculating the coefficients 
g^i, go, gi, ■ ■ ■ of (1.10). Our method is based on the expansion formula for the reflection 
coefficient, which was derived in [13]. This formula will be reviewed in section 3, after 
introducing some necessary notations in section 2. The calculation of gn is done in 
sections 4-9. 

It is easy to extend this method to potentials V{x) which have different asymptotic 
behaviors as a; — )■ — oo and x — )■ +oo. For example, wc can easily deal with the 
cases where V{x) approaches periodic functions with different periods as x ^ — oo 
and a; ^ +00. This will be discussed in section 10. 

If V{x) is asymptotically periodic, the corresponding Vs(a;) is also asymptotically 
periodic. That is to say, if V has the form of (1.7) with (1.8), then Vs has the form 

Vs{x) = V^{x) + Vl{x), V^{x + L) = V^ix), hm VHx) = 0. (1.12) 

But the converse is not necessarily true. Suppose that a Schrodinger potential Vs 
satisfying (1.12) is given. For simplicity, we assume that there arc no bound states. 
Then the corresponding Fokker-Planck potential V satisfies (1.7) and (1.8) only if the 
wave function ip{x) at the bottom of the lowest energy band remains finite for both 
X — > +00 and x — > — oo. This is what is called the 'exceptional case' in the conventional 
terminology of scattering theory [6]. In the 'generic case', the Fokker-Planck potential 
V{x) corresponding to an asymptotically periodic Vs{x) is not asymptotically periodic 
but tends to — oo at either x — )> +oo or x — >■ — oo (see example 2 in section 12). Our 
method is also applicable to Vs in the generic case, even though such Vs does not 
correspond to a Fokker-Planck potential satisfying (1.7) with (1.8). In the generic case, 
the expansion of Gs begins with the term of order (namely, g_i = in (1.10)). In 
section 11, we will see how to calculate go, gi, g2, . . . for the generic case. 
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Let the 2x2 matrix U{x,x'; k) be the solution of 
d 



dx 



U{x, x'] k) 



-ik f{x) 
f{x) ik 



U{x,x';k), U{x',x';k) 



1 
1 



We write the elements oi U as 



U{x, x'; k) 



a{x,x';k) I3{x,x';—k) 
f3{x,x';k) a{x,x';—k) 



(2.1) 



(2.2) 



and define the transmission coefficient r, the right refiection coefficient Rr and the left 
reflection coefficient Ri as 



t{x, x'; k) = 



1 



a{x, x'; k) 



a{x, x'; k) 



Ri{x,x'; k) = 



(3{x, x'; —k) 
q;(x, x'; k) 
(2.4) 



The generalized scattering coefficients r, Rr, Ri are defined with an additional variable 
W as 



f(x,x';W;k) 



^/l — C,'^ t{x, x'; k) 



Ri{x, x'; W; k) = Ri{x, x'; k) + 
where 

^{x, W) = tanh 



1 - ^Rr{x, X'- k) ' ^) - 1 _ ^R^^r,^ ^) ' 

iT'^{x,x'\k) 



1 — iRr{x, x'; k) ' 
W - V{x) 



(For an alternative definition, see (A. 2) and (A. 3) of appendix A.) We also define 

Rr{x,—oo;k) , , _ Ri{oo,x;k) 



Si{x,k) 



1 + Rr{x, — oo; k) ' 
S(x, k) = Sr(x, k) + Si(x, k). 

The Green function can be expressed in terms of this S as [14] 



1 + Ri{oo, x; k) ' 



Gs{x,y;k) 



1 



2ik^[l-Six,k)][l-S{y,k)] 
We use the notation 



exp 



ik{x — y) — ik / S{z,k)dz 



, 0-2, ■ ■ ■ , CTnll = ■■■ dzidz2 ■ ■ ■ dz„ cxp CrjV{Zj) 

J J a<zi<Z2<-<Zn<b '-j=0 



(2.5) 

(2.6) 

(2.7) 
(2.8) 

(2.9) 
(2.10) 



for n = 1,2,3,... and — oo < a < b < oo, where each aj is either +1 or —1. For 
simplicity, we write [+~]a) ^^'^ place of [+1]^, [+1, — l]a,! etc. In this paper, we 

also deal with integrals of the form of (2.10) with V replaced by Vp. We will denote 
such integrals with a left subscript 'p' as 



p[(7i,(72, . . . ,(7„]^ = j " j 



dzidz2 ■ ■ • dzn exp 



a<Zl<Z2<---<Zn<b 



■j=o 



(2.11) 
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As in [2] , we use the symbols 

1 . P 



P = A+]l-L. M^^_,, Lo^VpM, Vo^-log-, 

Q ^ A-+-+]:-L + p[+-+-]:-L- (2.12) 

All the quantities defined in (2.12) are independent of x. 
3. Formula for the expansion of Rr{x, — oo; W; k) 

In this section, we summarize the necessary results from [13]. We define the operators 
A and B, which act on functions of x and W, as 

Ah{x,W) = -^h{x,W), (3.1) 
ox 

d 

Bh{x, W) = -^ {smh[W - V{x)]h{x, W)} . (3.2) 
The generalized reflection coefficient Rr{x, —oo; W; k) satisfies the differential equation 

{A - 2ikB) [Rr{x, -oo; W; k) + ^(x, W)] - [l - ^^(x, W)] f{x). (3.3) 
Prom (3.3) we have 

Rr{x, -oo; W] k)^-i + {A- 2\kBy^{l - ^^)f, (3.4) 
where {A — 2ikB)^^ is the inverse of the operator A — 2\kB. This inverse is given by 

{A-2\kB)-'g{x,W)^ r dz ^ ^'^y^ ff(z,cD(x, z; W^; A;)), (3.5) 

J_oo 1- Rf{x,z;W;k) 

where 

-f w ^^ T/^ i + Ri{x,z]W;k) 

.{X, z; W; k) ^ V{z) + log , _ ^^^^^ ,y (3-6) 

We can formally expand {A — 2ikB)~^ in powers of k as 

{A - 2ikB)-^ = [1 + ik{2A-^B) + {ik)\2A-^Bf + ■■■ + {ikY{2A-^BY] ^"^ 

+ {ikf+\A - 2ikB)-^A{2A-^Bf+'A-\ (3.7) 
Substituting (3.7) into (3.4) yields the expansion 

^ = fo + i/cfi + (iA;)^f2 H h {ik)^fN + Pn (3.8) 

with 

fo^A-\l-e)f-i. (3.9a) 

f„ = (2^-^-B)"(fo + {n>l), (3.96) 

Pj, = (ifc)^+^(^ - 2ikB)-^AfN+i = 2{ik)''+\A - 2ikB)-^BfN. (3.10) 

Equation (3.8) makes sense if and only if the right-hand side of (3.96) makes sense for 
all n < N . (The remainder term p^v automatically makes sense if all f„ make sense.) 
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4. Inverse of the operator A 

In the formal expression (3.7), the symbol A~^ denotes the inverse of A. However, 
'inverse of A^ does not have a meaning unless we specify the domain of A . Specifying 
the domain of A amounts to specifying the boundary condition for Rr{x, — oo; W; k) at 
X — >■ — oo. Since the operator A in (3.3) acts on the function Rr{x, — oo; W] k) +^(x, W), 
the domain of A must be chosen so that -R^ + C belongs to that domain. 
We define the operators Aq^ and ^ by 

A-^^g{x,W)^ r giz,W)dz, (4.1) 

J —cx: 

-^p"^("'^)" Losinh(V-l-o) 

X J dz j d/|sinh[l/p(^') - W\g{z, W) - suih[Vp{z') - V^\g{z, Vo)]. 

(4.2) 

Let D/\ denote the set of two- variable functions h{x, W) which are picccwise 
continuously differentiable with respect to x, analytic with respect to W on the real 
axis and which satisfy 

lim h{x, W) = 0. (4.3) 

x-^—oo 

And let Dp denote the set of functions h{x, W) which are piecewise continuously 
differentiable with respect to x, analytic with respect to W on the real axis and which 
satisfy the conditions 

h{x + L,W)^h{x,W) and f Bph{z,W)dz^O for any x, (4.4) 

Jx-L 

where we have defined, corresponding to (3.2), 

Bph(x, W) = -^ {smh[W - Vp(x)]h(x, W)} . (4.5) 
It is easy to see that 

Ao^Ah = h if he Da. (4.6a) 
(We are allowing Ah to include delta functions.) And it was shown in [2] that 

Ap^Ah^h if he Dp. (4.6&) 

(In [2], Ap^, Bp etc are written without the subscript p.) In other words, if the domain 
of A is restricted to Da or Dp, the inverse of A is given by ^ or ^p^, respectively. 
If Rr{x, — oo; W; k)+^{x, W) belongs to Da, the domain of the operator A in (3.3) 

can be taken to be Da- Then the inverse of A is given by (4.1), and the expansion 
of Rr is obtained by letting A^^ = Aq^ in (3.9). This is the case for the non-periodic 
potentials discussed in [1]. On the other hand, if the potential is periodic, then Rr + ^ 
belongs to Dp (see [2] for details), and so A~^ should be used in place of A~^ in (3.9). 
In this paper, we are assuming that the potential has the form of (1.7). In order to 
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use the expansion formula shown in the previous section, we must find an appropriate 
domain of A for such asymptotically periodic potentials. 

Let denote Rr with V replaced by Vp. (That is to say, is defined in the same 
way as the definition of Rr in (2.1)-(2.5), with / in (2.1) replaced by /p.) We define 
R^ = Rr- RP, and 

U^,W)=t^nh U^,W) = ax,W)-CA^,W). (4.7) 

Then, we can write 

Rr+^ = R'r+^P + R^ + U- (4.8) 

Both R^{x, —oo;W; k) and ^a{x,W) vanish as a; ^ — oo, and hence it can be seen 
that R^ + e -Da- As for the periodic part, it was shown in [2] that R^ + G Dp. 
Therefore, Rr + ^ E Dp + D/s. (where h e Dp + D/^^ means that h — hp + hA with hp G Dp 
and /lA e -Da)- So we know that we should take Dp + Da as the domain of A. 

Prom now on, we assume that the operator A is defined with the domain Dp + Da- 
To calculate (3.9), we need the inverse of A. li h E Dp + Da and g — Ah, we can write 

9^9p + 9A, (4.9) 

where gp — Ahp and gA — Ah a with hp G Dp and /ia £ -Da- Any function g belonging 
to the range of A can be uniquely decomposed into two parts as (4.9), where gp is a 
function satisfying 

gp{x + L,W)=gp{x,W), f gp{z,W)dz = for any x, (4.10) 

Jx-L 

and gA is a function such that 

/X 
5'a(-2, W) dz exists and is finite. (4-11) 
-oo 

It is obvious that gp satisfies (4.10) if gp = Ahp with hp G -Dp. Conversely, if gp 
satisfies (4.10), it can be shown that A~^gp G .Dp and AA~^gp = gp (see [2]). From 
A'^g = hp + hA and (4.6), we have 

A-'g = Ap^gp + A^^gA. (4-12) 

Thus, when the potential is asymptotically periodic, the expansion of Rr is given by 
(3.8) and (3.9) with A^^ acting as (4.12). To calculate A~^g, we first express g as the 
sum of the periodic part gp and the non-periodic part gA, satisfying (4.10) and (4.11), 
respectively. Then, the right-hand side of (4.12) is calculated with the operators defined 
by (4.1) and (4.2). 

5. Expressions for f„ 

Now let us calculate the coefficients f„ given by (3.9). To calculate fg, it is necessary to 
decompose (1 — ^^)/ into periodic and non-periodic parts. Note that 

[1 - wm.) = —zi^iLyM = w). (5.1) 
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Using (4.7), we write (5.1) as 

[1 - W)]f{x) = ^U^, W) + ^U{x, W). (5.2) 

The two terms on the right-hand side of (5.2) correspond, respectively, to Qp and Qa of 
(4.9). They satisfy conditions (4.10) and (4.11). So, according to (4.12), 

^-'(1 - e)f = ^p'^U^^ + A^'-^U^, W). (5.3) 
The first term on the right-hand side has aheady been calculated in [2]. The result is 

A'-^U^^ = ^p(^' ^) - t^nh (5.4) 
(see equation (8.1) of [2]). It is obvious that the second term of (5.3) is 

A'o^Uix, W)^ j dz —U{z: W) = W) - e(^, W) - ep(x, W). (5.5) 

(Note that A^JI^a = but ^p^^^p ^ since e -Da but ^p ^ Dp.) Therefore, 

^-'(1 - e)f = - tanh "^^^ + e(x, W), (5.6) 
and (3.9a) gives 

W -Vq 

ro = - tanh ^ . (5.7) 

Let us proceed to the calculation of fi = 2^~^B(fo + C)- can be shown that 

(see equation (8.3) of [2]). We decompose the right-hand side as 

2B{fo + 0=<ll + <lt (5-9) 

where 

^0 = cooh^^^-^" sinh[K) - Vp{x)l = ^^^J'w-Vp sinh[^ - V{x)] - q^. (5.10) 

It is easy to check that satisfies conditions (4.10). In order that q^ satisfy (4.11), 
it is necessary that V{x) approach Vp{x) sufficiently rapidly as a; — >■ — oo (see the next 
section). Assuming that this condition is satisfied, we have, from (5.9) and (4.12), 

2A-'B{fo + = + A'Qo- (5-11) 

So we obtain fi as a sum of two terms 

fi=fP + ff, (5.12) 

where 

1 

cosh^M 

A' ^^^J^_yJ smh[Vo - V{x)] - sinh[yo - Vp{x)]}. (5.136) 



^i=A^ — r,2w-Vo smh[Vo - Vp{x)], (5.13a) 
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The right-hand side of (5.13a) can be calculated using (4.2). Details for the calculation 
is given in [2]. As a result, we have 

(See equation (8.7) of [2]. Note that [ai, . . . , an\a in [2] is p[(Ji, ■ ■ ■ , cr„]^ in this paper.) 
Equation (5.136) can be written as 

r-f = ^^^Jw^ - e-^oA+(.)] d., (5.15) 

where we have defined 

A±(x) = e±^(^) -e^^''^^). (5.16) 

Obviously vanishes as x — > — oo (provided that the integral on the right-hand side of 
(5.15) is convergent), while is a periodic function of x. The first-order coefficient fi 
is thus obtained as the sum of the periodic part (5.14) and the non-periodic part (5.15). 

To calculate f„ for larger n, we can use the recursion relation f„ = 2A~^Bfn-i-i 
which follows from (3.96). We define 

^„ = 2Bf„ (n>l), (5.17) 

so that 

f„ = A^^qn-i. (5.18) 
We assume that f„ and g„ can be written as the sum of periodic and non-periodic parts, 
rn^fl + ft qn^ql + qt (5-19) 



where 



fl{x + L]W) =fl{x,W), lim f^(a;,iy) = 0, (5.20a) 

ql{x + L-W) = ql{x,W), lim q^{x,W) = Q. (5.206) 

This assumption will be justified by the result. We split B into two parts as 

B = Bp + Ba: (5.21) 

where Bp is defined by (4.5), and Ba{= B — Bp) can be expressed in terms of A='= 
(equation (5.16)) as 

e^ / a \ e"^ / d \ 

= ^A-(.) (l + -) + — A+(x) (l - ^) . (5.22) 

From (5.17), (5.19) and (5.21), it follows that 

ql = 2Bp fP, q^ = 2Bf^ + 2Ba (5.23) 

It can be shown [2] that q^ satisfies conditions (4.10), and so Ap^q^ makes sense. 
Assuming that A^^q^ also makes sense, we have A~^qn = ^p^Qn + ^o^Qn- Therefore, 

rl = 'C-i = fP_i, = A^'qti = 2A^'B fti + 2 A"'Sa (5.24) 
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By iterating (5.24), we obtain 

= (2A;'B,r-'fl (5.25) 

n-2 

rt = {2A^'Br-'f^ + J2i2A^'By 2A^'Ba fl_,_,. (5.26) 

i=o 

In this way, and can be calculated from and ff (equations (5.14) and (5.15)). 
It is obvious from (5.25) and (5.26) that equations (5.20a) are satisfied as long as the 
right-hand side of (5.26) makes sense. If Va = 0, then f„ = f^. The are identical 
with the f„ studied in [2] for periodic potentials. 

Let us explicitly write oTit the expression for f2 = + f^. The calculation of f\ 
was done in [2]. The result is given by equation (8.18) of [2] as 

'^^"4Locosh^M\^ p[+ e p[ + J._^ 

.— (Lj + 4Q)sinh — 

with Q defined by (2.12). For n = 2, the right-hand side of (5.26) consists of the 
two terms 2A^^Bf^ and 2A^^B^fl. Applying (5.22) and (3.2) to (5.14) and (5.15) 
respectively, we can calculate 

e(w^+Vb)/2 A-fxl +e-(^+^o)/2A+('x) / \ 

Bft = 4eosh^M — J ^ [^''"^"(^) - (5-286) 

Hence, the non-periodic part is obtained as 

1 



f — {Lt + 4g) sinh -—^ , (5.27) 



^2 



A 



4LoCOsh3i5^ 

X /_'^{[e(^-^^°)/^A-(.)+e-(^+WA+(.)] {a+'Y.-l - A'+Y.-l) 

+2Lo(e(^+W[_]^ + e-(^+^°)/2[+]^) [e^°A-(^) - e-^°A+(^)] | d^, 

(5.29) 

where we have used dz dz' = dz' dz. 



6. Condition for the finiteness of r 



n 



Equation (3.8) is meaningless unless the coefficients f„ = r^-l-f^ are finite for all n < N . 
The periodic part f^, which is given by (5.25), is finite for any n (see [2]). However, the 
terms on the right-hand side of (5.26) are not necessarily finite. For these terms to be 
finite, it is necessary that V/^{x) vanish sufficiently fast as a; — > — oo. 

Since A^(z) ~ ±e^p'^^)VA(-2) as 2; — > —00, and since V^{z) is a bounded periodic 
function, there exists a constant C such that 

|A±(^)| <C|Va(^)| (6.1) 
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for — oo < z < X with fixed x. Using this in (5.15), we find 

|rf I < C I \V^{z)\dz. (6.2) 

J — oo 

(We will use the symbol C to denote a finite constant which is not necessarily the same 
at each occurrence. We regard x and W as fixed.) We can see from (5.286) that 

|2^o'^^f I <C I dz f dz' |e^"A-(z') - e-^°A+(z')| , (6-3) 

J —oo J —oo 

since e^^*^^) are bounded for — oo < z < x. Using (6.1) in (6.3) gives 

/X PZ PX PX 

dz dz' \Va{z')\ = C dz' dz\VA{z')\ 
■oo J —oo J —oo J z' 

r (x-z)\VA(z)\dz. (6.4) 

J —oo 

More generally, it is easy to see that 

|(2A-^B)"-^ff I < C|(A-')"-'rf I < CiAo'r-' f \VA{z)\dz. (6.5) 

J —oo 

Hence we can derive, in the same way as (6.4), 

pA^^BY-^f^\ <C r (x- zY~^ \VA{z)\dz. (6.6) 

./ — oo 

This gives an upper bound for the first term on the right-hand side of (5.26). Similar 
inequalities hold for the remaining terms of (5.26). Using (5.22) and (6.1), and also 
using the fact that |fJJ_j_i| < C and \-^f'^n_j_^ < C, we have 

\2Ao'BAfl_j_,\<C r \VA(z)\dz, (6.7) 

^— oo 

and, just like (6.6), 

\{2A^^By2A^^BAfl_j_,\<C r {x-zy\VA{z)\dz. (6.8) 

J — oo 

Let F^"-* denote the set of functions v{x) which satisfy the condition 

(1 + |x|"')|v(x)| dx < oo for any finite a. (6.9) 

The right-hand side of (6.6) is finite if Va G -^i-i- Then, the right-hand side of (6.8) is 
also finite for j < n — 1. So, we know from (5.26) that \f^\ < oo if Va € -^ili- Since 
is always finite, it follows that 

\fn\ < OO if Fa e F^ZI (6.10) 

If |fiv| < oo, then |f„| < oo for all n < N. Therefore, we can conclude that (3.8) makes 
sense if Va e -^jv-i- 



J — c 
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7. Small-fe behavior of the remainder term 

In the previous section, it was shown that equation (3.8) makes sense if Va G FnIi- In 
this section, we will show that = o{k^) as A; ^ if Va € p'n-i- 
Using (3.5) and (5.17), we can write (3.10) as 

Pn = {ik)''^^ j Q{x, z; W; k) qN{x, 6j(x, z; W; k)) dz, (7.1) 

J — oc 

where 

l-[Ri(x,z;W;k)r ^^'^^ 
Using (5.19), the right-hand side of (7.1) is spht into two parts as 

Pn = {ikf+^ j Q{x, z- W- k) (g{z, u{x, z- W- k)) dz 

+(iA;)^+^ /" Q{x,z;W;k)q^{z,u{x,z;W;k))dz. (7.3) 
Let rP and Rf denote, respectively, f and Ri with V replaced by Vp. We define 



\f'P{x,z;W;k)f 



(7.5) 



In the limit A; — > 0, we have the expressions 

W_V(z) _ W-VJz) 
lim t(x, z; W; k) = sech —, lim tP(x, z: W; k) = sech 

k-^O 2 k-^O 2 

lim Ri{x, z; W; k) = tanh KjzXl^^ Um RP(x, z; W; k) = tanh ^ ~ ^^^^^ (7.6) 

k — ^0 2 k — ^0 2 

(see equations (2.18) of [2]). Hence, 

lim Q{x, z; W; k) = lim QJx. z; W; k) = 1. (7.7a) 

lim uix, z; W; k) = lim uJx, z\ W; k) = W. (7.76) 

It can also be shown that 

lim / Q{x,z;W;k)q^j^{z,uj{x,z\W\k))dz 

/X 
Qp{x, z\ W; k) f^{z, cdp{x, z; W; k)) dz, (7.8) 
-oo 

provided that Va G Fq (See appendix A for the derivation.) It was shown in [2] that 
lim r Qp{x,z;W;k)q%iz,Upix,z;W;k))dz ^ A;'f^{x,W). (7.9) 

The right-hand side of (7.9) is equal to f^^i (see (5.24)). Therefore, 

lim / Q{x,z;W-k)fj,{x,u{x,z-,W-k))dz^f''j,,^{x,W). (7.10) 

J-oo 
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This gives the small-/c behavior of the first term of (7.3). 

Next, we consider the second term of (7.3). Substituting (5.26) into (5.23), we write 

iV-l 

QN^U+^Vm + W, (7.11) 
m=l 

u = 2B{2Ao'Bf''ft vm = 2B{2Ao'l3r-' 2A^'Ba fl_^, w = 2Ba f%. (7.12) 

As shown in {7.7b), the quantity u}{x, z; W; k) tends to 1^ as A: — > 0. This approach is 
uniform in z. Wc fix x and W, and let be a fixed (sufficiently small) positive number. 
Using the same argument as in the last section, wc can easily show that 

\u{z,(:j{x,z;W;k))\ < ChN{z), \vm{z,Lu{x, z;W; k))\ < Chm{z), 

\w{z, u{x, z]W]k))\< Cho{z) (7.13) 

ior —oo < z < X and \k\ < ko, where we have defined 

hr^{z)= r {z-zT-'\VA{z')\dz' (n>l), ho{z) = \Va{z)\ . (7.14) 

J —oo 

It is easy to see that hn{z) dz < oo if Va e Fi~\ We have the inequality for Q, 

\Q{x,z;W;k)\<l, (7.15) 
as shown in appendix B. From (7.11), (7.13) and (7.15), we find that 

N 

\Q{x,z;W;k)q^{z,u{x,z;W;k))\<h{z), h{z) = C hm{^) ■ (7.16) 

m=0 

Obviously, J^^h(z) dz < oo if Va £ Fn~^- Therefore, by the dominated convergence 
theorem, the integral in the second term of (7.3) commutes with the limit A; — >■ if 
Va e Namely, if Va G Fj^^ we obtain, by using (7.7), 

/X 
Q(x, z; W; k) q^(x, uix, z; W; k)) dz 
-oo 

lim Q(x, z; W; k) q^(x, u}(x, z; W; k)) dz 
f q^{x, W) dz = Ao'q^ix, W) = f^^,{x, W). (7.17) 

J — oo 



Prom (7.3), (7.10) and (7.17), we have 



lim jr^jWT^PN = r^N+i + ^n+i = ^n+i if I^a e Fj^, \ 



(7.18) 



Since — {ik)'^'^^rN+i + Pn+i, equation (7.18) implies that limjfc^oPAr+i/(i^)'^''"''^ = 
if Va e Replacing by - 1, we obtain 

lTo(i^^^ = ° Vke41. (7.19) 
Therefore, if Va e F^^}^, then can be expanded to order /c^ as 

Rr{x, -oo; W; k) = fo + ikn + (iA;)¥2 + ■ ■ ■ + (i/c)^fjv + o{k^). (7.20) 
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8. Expansion of S 

It is shown in [1] that the expansion of Sr (defined by (2.7)) is obtained from (7.20) as 
Sr{x, k)-^^ a^{x) + ika^ix) + {ikyaf{x) + ■■■ + {ikfa^{x) + o{k^), (8.1) 
where 

a^{x) = - hm e-^+^(^) \f,{x) - 1] , (8.2a) 

a^{x) = \ hm e-^+^(^¥„(x) (n > 1). (8.26) 

Equation (8.1) is vahd if Va £ -^jvl\, as is (7.20). Putting into (8.2) the expressions for 
fo, fi and f2 (equations (5.7), (5.14), (5.15), (5.27), (5.29)), we obtain 

= -^e^(^)-^°, (8.3a) 

2 J-OO 

= ^e^(^)-^o {e-^%[+-+]:_. - ^ W + 4g) } 

j^^v{x)-2v, r dz{e^°A-{z)-e-^°A+{z))[+Y^. (8.3c) 

^ — oo 

The expansion of Si can be obtained in the same way. Let F^^"* denote the set of 
functions v{x) satisfying 

/"OO 

/ (1 + |a;|"')|t'(a;)| dx < oo for any finite a. (8.4) 

J a 

If Va e Fi^}^, we have 

5;(a;, k)-- = a^{x) + iA;a5^(a;) + {ikfa^{x) + ■■■ + {ikfa)^{x) + o(A;^), (8.5) 
2 



where 



«o(^) = -ie^(^)-^°, (8.6a) 
2 Jx 

"21^"''^'^''''" ^^^^^ (p[+-]U - p[-+]^-l) 

/■oo 

^gy(x)-2i-o / (e^o^-(^) _ e-^oA+(z)) [+]^. (8.6c) 
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If Va G F^-\ n F^+i, i.e. if Va G L],- wc can add together equations (8.1) and 
(8.5) to obtain the expansion of S (equation (2.8)) as 

S{x, k)-l^ So{x) + iksi{x) + iik)^S2{x) + ■■■ + (iA;)^s^(x) + o{k^) (8.7) 
with 

Sn{x) =a^{x) +ali{x). (8.8) 
From (8.3) and (8.6) we have 

so{x) = -e^(-)-^o, (8.9a) 

„V{x)-Vo poo 

s^{x) = —— / dz PA-(z) - e-^°A+(z)] , (8.96) 

J —oo 

POO 

- dzU+-]U-,[-+]U)A+{z) 

J X 

+2Lo r d;^ [e^«A-(;2) -e-^«A+(^)] [+]^ 

+2Lo^ d^ [e^oA-(^) -e-^°A+(z)] [+]^|. (8.9c) 
In a similar way (we omit the calculation) , we can derive the expression for S3 as 

V{x)-2Vo ( roo / r3 0O\ 

+ r dz{,[+-]U -,[-+]:_,) [3e-^<'A+(z)-e^''A-(^)] [+]^ 
J —00 

roo 

- / dz U+-]U - [3e-^°A+(^) - e^°A-(^)] 

J X 

+2Lo r d^ [e^°A-(^) -e-^°A+(^)] (3e-^°[++]^ - e^°[-+]^) 
J —00 

+2Lo [ dz [e^"A-(;2) - e-^°A+(^)] (3e-^''[++]^ - e^°[+-]^) j. (8.9 (i) 



9. Expansion of the Green function 

The expansion of the Green function in terms of k can be obtained by substituting 
the power-series expression of S into (2.9). To expand Gs to order k^ , we need the 
expansion of -S" to order k^'^^ (as shown below). Let us define 



px 

(ln{x,y) = - / Sn-l{z)dz, tn{x) = 



(•^) 

so{x) 



-e^»^^(")s„(a;). (9.1) 
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If Va e Ljy, then equation (8.7) holds with N replaced by A'" + 1. Substituting this into 
(2.9), we obtain equation (1.10) with 

g-i{x,y) = \ (9.2a) 

2^^ SQ[x)so[y) 



9o{x,y) = 



ti{x) ti{y) 
QiM--^ ^ 



9~i{x,y), {9.2b) 



gi{x,y) = U2{x,y) + — [tx{x) + ti{y)\ 



2 2 8 8 4 



^2(2;,?/) = ig3(2;,2/) + gi(x,y)g2(2;,z/) + 



6 



M^ + ^)[t,(.)+t,(,)] 



^) / 3[ti(x)]^ ^ 3[^i(y)]^ ^ ^i(^)^i(y) ^2(x) hjy) 



tsix) t^jy) ^ StijxMx) ^ 3ti{y)t2{y) 5[ti{x)f 5[ti{y)f 



2 2 4 4 16 16 

and so on. It is easy to see that we need So,Si, . . . , sjv+i in order to calculate g^. Thus, 
(1.10) holds if Va £ L]^- Prom (8.9) and (9.2), the explicit forms of g^i and go are 
obtained as 

^_l(x,|/) = ^e-[^(^H^(^)l/^ (9.3a) 



^o(x, y) = i I [+]^- +"^l^dz [e^o A-(z) - e"^" A+(z)] | e-[^(-)+^(^)l/^ 



(9.3&) 



10. More general potentials 



The method presented in this paper is also applicable to the cases where V{x) does 
not approach the same periodic function as x — )■ —00 and x — )■ +00. Suppose that 
V{x) Vpi{x) as a; — )■ —00 and V{x) — )■ Vp2{x) as x — > +00, where Vpi and Vp2 are two 
different periodic functions with periods Li and L2 respectively. Namely, 

V = Vp, + Vai = Vp2 + Va2, Vp,{x + Li) = Vp,{x), Vp2{x + L^) = Vp2{x), 

lim Vai{x) = 0, lim Va2{x) = 0. (10.1) 

The expressions for and given in section 8 (equations (8.3) and (8.6)) still hold 
if we use Vpi for and Vp2 for a^. The expansion of Gs is obtained by calculating 
Sn — + a^, and substituting into (9.1) and (9.2). 

Let p^o"!, . . . , 0"rJa, Vot, i^oi ^^d (i = 1, 2) be the quantities defined in the same 
way as p[cri, . . . , cr„]^, Vq, Lq and A''^, respectively, with Vp, Va and L replaced by Vpi, 
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VAi and Li. If Vai G f}^^ and Va2 G F^^l, the Gr een function can be expanded to order 
in the form of (1.10). Instead of (9.3), we have 



where 
ii(x) = — 



-1 



2L 



02 



(p2[+-]^-L.-p2[-+]^-L.)+e-^- r d^[e^-Ar(^)-e-^-A+(^)] 

+e-^«2 ^ d;2 [e^"^A2 (;2) - e-^"^A+(;2)] |. (10.3) 

It is equally easy to calculate the expansion of the Green function for the cases in 
which V{x) is asymptotically periodic as x — )■ — oo and V{x) ±oo as x — >■ +oo. For 
such cases, we can use (8.3) together with the expressions for given in [1]. 

11. Generic case for the Schrodinger equation 

Let us suppose that the Schrodinger equation is given with an asymptotically periodic Vs 
satisfying (1.12), and that the corresponding Fokker-Planck potential V is yet unknown. 
Wc assume that there are no bound states, and we shift the origin of the energy scale 
(/c = 0) to the bottom of the lowest band. The Schrodinger equation at A; = reads 

-^V'o(x) + V^{x)i,o{x) = 0. (11.1) 

Let iPq{x) and iPq{x) be the solutions of (11.1) which remain bounded as a: — )■ +oo 
and X — 7- — oo, respectively. We can choose the phase of so that V'o^(^) > for 
any x. (Although if)"^ still has an arbitrariness of a positive multiplication factor, this 
arbitrariness does not remain in the final result for the Green function.) Let us define 

V^{x) ^ -2\og^l^^{xl f^{x) ^ -\^V^{x). (11.2) 

It is easy to check that both and satisfy (1.4). So both V'^ and V~ are 
Fokker-Planck potentials corresponding to V^. Since V^{x) and V~{x) tend to periodic 
functions as x — >■ +00 and x — >■ — 00, respectively, we can write — + V^, where 

V^{x + L) = V^{x) \imV^{x)^Q, lim VX(x) = 0. (11.3) 



If ipQ and iIjq coincide up to a multiplicative factor, then = V~ + C with a 
constant C, and so V^{x) are asymptotically periodic for both x — )> +00 and x — )■ —00 
(that is, we have limx^_^V^{x) = and lima;_j.+oo (a;) = in addition to (11.3)). 
This is the so-called exceptional case. In this exceptional case, we can use the result of 
section 9 for the Green function, using either y+ or V~ in place of V. On the other hand. 
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if ipo and are linearly independent, then V~^{x) and V~{x) arc not asymptotically 
periodic for x — )■ — oo and x — )■ +00, respectively. This is the 'generic case'. 

To deal with the generic case, we can use the method described in section 7 of [1]. 
As explained there, the expansion of S takes the form 

S{x, k) -1 = (iA;)"^s_i(a;) + so{x) + iksi{x) + {ikfs2{x) 

+ --- + {ikfsN{x) + o{k'') (11.4) 

with 



S-l^-[f-(x)-f^x)], Sr, 



{n > 0), 



(11.5) 



where a^^ and 0^=*= denote, respectively, and with V replaced by V^. Equation 
(11.4) is valid if e and e Fj^\. This is a generahzation of (8.7). In the 

exceptional case (where = V~ + C), we have a^~^ = a^~, aj^+ = ajj~ and s_i = 0, 
and hence (11.4) reduces to (8.7). If ipQ are hnearly independent, then f~^{x) ^ f~{x) 
for any x, as can be easily verified. So, in the generic case, s^i{x) ^ for any x. 

In the generic case, the expansion of Gs to order k'^ is obtained from the expansion 
of S to order k^^^. Substituting (11.4) (with replaced by — 1) into (2.9), we have 

Gs{x, y; k) = go{x, y) + \kg^^{x, y) + (\kfg2{x, y) H h {ik)^gN{x, y) + o{k^), (11.6) 

where 



g2{x,y) 



go{x: y) = 

(l2{x,y) + 



-exp[go(a:,y)] 
2v/^s-i(x)s_i(2/)' 

Qiyx^y) 



2s-i(a;; 
[(li{x,y)f qi{x,y) 



+ 



8 



so[x) 
s-i{x) 



2 

1 2 



_ so{y) 

2s-i(i/) 
So(a;) 



^0(2;, t/), 



+ 



+ 



8 



sojy) 



+ 



s_i(x) ' s_i(y) 
so{x)so{y) 
4s_i(x)s_i(y) 



si(a;) 



+ 



(11.7a) 
(11.76) 

siiy) 



and so on, with Qn defined by (9.1). For A" > 2, equation (11.6) is vahd if e 



s_i(x) s_i(y) 

9o{x,y), (11.7c) 

{-) 



and y+ e 



The explicit forms of go and gi obtained from (11.7), (11.5), (8.3) and (8.6) are 
- exp {[V-{x) - V+{x) - V-{y) + V+{y)] /4} 



9i{x,y) = 2 



- /n^)] [f-(y)-fHy)] 

qV-{x)-V- ^V+(x)-V+ QV-iy)-V- _^ ^V+{y)-Vo 



11.8a) 



+ 



^o(a;,y). 



;il.86) 



f-{x)-f+{x) f-(y)-f+(y) 

where V^q^ and Vq~ are defined by (2.12) with Vp replaced by and V~. respectively. 
By substituting (11.2) into (11.8a), we can reproduce the well-known expression for gQ-. 

'^oix)i^o{y) 



go{x,y) 

where W denotes the Wronskian defined by W[ijj, 0] = — ip'cj). 



(11.9) 
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Figure 1. The potential V{x) of example 1 (equations (12.1)). 



12. Examples 

Let us demonstrate the calculation of the expansion with simple examples. 
Example 1. 

We consider the potential V{x) = Vp{x) + Vi\{x) with 

(0 < X < a) 

/ Vp{x + L) = Vpix), (12.1a) 

C [a < X < L) 




-h (0 < X < a) 
(otherwise) 



[12.1b) 



This potential is illustrated in figure 1. The Green function can be exactly obtained for 
this potential. Before showing the calculation of the low-energy expansion, let us first 
explain about the exact Green function and its properties. 

For any potential, the Green function can be expressed in terms of r, Rr and Ri as 

Q^^^y.f.^^ [1 + ^;(oo, x; k)] [1 + Rrjy, -oo; k)] t{x, y; k) 2) 

' ' 2ik\l — Ri{(yD^x]k)R^{x,—oo]k)\[l — Ri{x,y]k)Rr{ii,—oo]k)] 

(see equation (3.6a) of [14]). We confine ourselves to the case Q < y < x < a. Then, for 
this potential, Ri{x,y; k) = and T{x,y; k) = e^^^^^y\ So, (12.2) becomes 

Gs(x, y; k) = + ^li^^ k)] [1 + R^jy, -00; k)] e'^(-^) ^ ^^^.3) 
' 2ik[l — Ri {00, x; k)Rr{x, —oc] k)] 

The exact expressions of the reflection coefficients are 

Rr{x, -00; A:) = ^2i/c. ^tanh(V^2j^^^^^^ ^^^^^ ^_ ^ ^^^^ _ ^^^.4) 

for < X < a, where 

aik) ^ (1 - Ah'^'^") , m = {^,A (e^^^^ - l) , tanh -f, 

_ -a{k) + a{-k) - iv/4 - [a{k) + a{-k)]^ 
^° = WFk) ■ ^ 
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Figure 2. The real and imaginary parts of Gs{x,y;k) for tlie potential shown in 
figure 1, plotted as functions of real k, with various values of h. (a) h = 0, (b) h = 0.2, 
(c) h = 0.5, (d) ft, = 1. In all the graphs, C = 1, L = 1, a = 0.6, x = 0.4 and y = 0.1. 
The graphs are shown here for the range of k in the lowest two energy bands and 
the lowest two gaps. The bands are < k < ki and k2 < k < k^, and the gaps are 
ki < k < k2 and k^ < k < k^, where fci ~ 2.21, fc2 ~ 4.02, k^ ~ 5.77 and k^ ~ 6.88. 
The imaginary part of Gs is identically zero in the gaps. 



(We omit the derivation of these expressions.) The exact Gs is given by (12.3) with 
(12.4). The graphs of the exact ReGs and ImGs, as functions of real k with fixed x 
and are shown in figures 2 and 3. The case h = corresponds to a purely periodic 
potential, which was studied section 11 of [2]. When h = 0, the graphs are discontinuous 
at the edges of the energy bands (figure 2(a)). The graphs become continuous when 
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Figure 3. Same as figure 2, with larger values of h. (a) h — 2, (h) h = 4, (c) h = oo. 
In (b), the real part of Gs(fc) tends to a finite value (~ 1000) as fc ^ 0, although 
the graph is truncated at the top. The graph in (c) is the plot of (12.6). This graph 
diverges to infinity as — >■ 0. The singularity at k = vr/a ~ 5.24 corresponds to an 
eigenvalue for the infinite square well potential. The peak around k = ir/a in the 
graphs of ImGs becomes a delta function as we let h — >■ +00. 



h ^ 0, but singularities remain at the edges of the bands (figures 2 (b)-(d)). As h is 
further increased, these singularities become less prominent in the graphs (figures 3(a) 
and (b)). In the extreme case h = +00, the potential is an infinite square well, for which 
the exact Green function is 

Gsix,y;k) = '''^^^l:f'^'^^y^ (12.6) 

(see figure 3(c)). As h becomes large, the Green function approaches (12.6) except at 
k = rnr/a {n integer), where the right-hand side (12.6) has poles. 

Now let us turn to the low-energy expansion. Since Va G L]^ for any N, expansion 
(1.10) is valid for any A^. We need to calculate the various quantities appearing in (8.9). 
The quantities involving only Vp (and not Va) have already been calculated in [2]. It is 
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shown in section 11 of [2] that 

Lo = V(« + &e^)(« + &e-^), yo = \ log "^/'^g , 

2 a + oe~'" 
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(12.7a) 
(12.76) 



p[+-]"-L = i^{o.' + 6") + 6-^^6(0 - x) + e'^fex, 



(12.8a) 



p[+-+]x-L = + 3^ ) + - «) sinhC + -(3a' + . (12.86) 

The expressions for p[ — and p[ — I — obtained from (12.8a) and (12.86), 
respectively, by changing the sign of C. Other integrals in (8.9) are also easy to calculate. 
For < z < X < a, we can see that 

„-2/i Y 



[+]: = e-\x-z), [++] 



- [-+T. = [+-] 



{x-zy. (12.9) 



The functions A''^ defined by (5.16) are 

A-WJ«''-1 (0<.<a) 
I (otherwise) 

The integrals including A='= can be calculated, for example, as 

/oo 
dz [e^°A-{z) - e-^°A+(^)] = 2a [sinh(Vb + h) - sinh Vq] 
-oo 

Substituting the above expressions in (8.9), we obtain 

So = -e-^^°+''\ si = e-^^°+^^a[smh{Vo + h) - sinh V"o], 



;i2.10) 



:i2.111 



S2 = 



-(a^ + 36') + ^(3a' + b^)e'^ + 2bx{x - a)e-^ sinh C - ^ (L^ + 4Q) 
6 6 8ivo 



+^-2(Vo+h) [^2 ^ 

g-(2yo+ft) 



a — 



^] [sinh(K) + ^) -sinhT/o], 



2Ln 



(3e-vb + e^°) ^(a^ + 36^) + (3e-^°+^ + e^""^) ^{3a' + 



-ba^ (e-'^ + l) e-^«sinhC 



1 



2Lr 



+ 



Q~i2Vo+h) 



(X — a 



Vo-2h Vo 



e^«) [sinh(yo + h) - sinh Vq] , 



^12.12) 



with Lq, Vq and Q given by (12.7). (We have used Lq = 6sinhC/ sinhVo-) Substituting 
(12.12) into (9.1), and then into (9.2), we obtain g_i, qq, gi and g2- In particular, 

QVo+h I 



9-1 = 



9o = ^[x-y) + ^ e^°+'* [sinh(\/o + h) - sinh , 



;i2.13) 



as can also be seen directly from (9.3). Thus, the expansion of the Green function is 
obtained to order k'^. This result is shown in figure 4 (the broken lines). We can see 
from these graphs that this expansion is indeed correct. 
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('^) ReGs ImGs 




Figure 4. KeGs{x,y;k) and lmGs{x,y;k) plotted as functions of real k, for (a) 
h = 0.1, (b) h — 0.5 and (c) h — 1. (All the other parameters are the same as in 
figures 2 and 3.) The solid lines show the exact Green function (close-up of figure 2). 
The broken lines show the result of the low-energy expansion up to order k^. The 
dotted lines represent the approximation given by (12.21). 



For a potential like (12.1), there is a good approximation method which can be 
used for a wide range of k. Let us explain this approximation in connection with the 
power-series expansion of Rr- From (5.7), (5.14), (5.15), (5.27) and (5.29), we have 



ro = tanh 



Vo-W 



cosh^ ^ 



X sinh(Vo + h) — - sinh Vq 



:i2.14a) 



r2 



cosh' 



ppjT— [x^ sinh(Vo + h) — ax sinh Vb] cosh 

9. 



W + Vo + 2h 



l2Lr 



(a^ + 362) g-^j^ 



W + Vn 



12L, 



-{3a^ + b^)smh 



W + Vo-2C 
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+1^ + ^^^^ ^^-y^}- (12.146) 

Recall that Rr{x, — oo; A;) is obtained from Rr{x, —oo; W; k) by setting W = V{x). Since 
V{x) — —h for < x < a, the expansion of Rr{x, — oo; /c) for < x < a is 

Rr(x, -oo; A;) = ro + i^n + (ifc)V2 H , (12.15) 

where ro and ri are obtained by letting W — —h in (12.14a) as 

To = tanh ° ^ — , ri = 2(a; — 5)ro, (12.16) 
with 5 defined by 

asinh\/o 

'^=2sinh(yo + /.)- ^''-'^^ 
Setting = — /i in (12.146), we can see that 

r2 = 2(x - 5) Vo + • • • , (12.18) 

where the terms represented by the dots behave like e"'*/^, and hence are negligible, 
when h is large. In the same way, it can be shown that 

on 

r„~-(a;-5)"ro (12.19) 
n\ 

for any n, when h is large. Substituting (12.19) into (12.15) gives 

Rr{x, -oo; k) ~ roe^^'=(^-'^). (12.20) 
Prom (12.3), (12.20) and the second equation of (12.4), we obtain the approximation 

In fact, this approximation is equivalent to replacing the potential V{z) by the effective 
square well potential 

, , [ {6 < z <a-6) 

VMz) = ^ (12.22) 

I Vq (otherwise) 

The result of this approximation is also plotted in figure 4 (the dotted lines). As we can 
sec from figures 4(b) and 4(c), equation (12.21) gives a good approximation, when h is 
not very small, for a wide range of k except near the band edges. This approximation 
becomes better for larger h, so that the graphs shown in figures 3(a) and 3(b) can be 
very accurately approximated by (12.21). 

Example 2. 

Next, we study an example where the Schrodinger potential Vs, rather than the Fokker- 
Planck potential V, is given. We consider Vs{x) = Vp{x) + V^{x) with 

„ f —Eq (0<x<a) „ Q 

y^i^) = {^ p V^\x + L) = V^{x), {12.23a) 

[C - Eo {a<x < L) 
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Vs{x) 




Figure 5. The Schrodinger potential Vs of example 2 (equations (12.23)) and the 
corresponding Fokker-Planck potentials V^. The parameters used for these graphs 
are C = 1, i = 1, a = 0.6 and h = 0.5. From (12.24) we have Eq ~ 0.3952. 



V^{x) = { ^ . (12.236) 

(otherwise) 

(See figure 5.) The periodic part is a Kronig- Penny potential [15]. We assume /i > 
so that there are no bound states. The constant Eq is determined by the condition that 
the energy at the bottom of the lowest band be zero. It is easy to see that Eq is the 
smallest solution of the equation 



V^C^tanh b[£E^ = /^tan (12.24) 

The expansion of the Green function can be obtained by using the method explained in 
section 11. We define 

q .bqp ap 



p=^Q, q=^C- Eo, s^^h- Eo, ^ = ^tanh^ = ^^tan^. (12.25) 

s 2 s 2 

The solution of (11.1) which remains bounded as x — ?■ +oo is 

{cosh [s{x — a)] — ^ sinh [s{x — a)] (0 < x < a) 
sechf cosh[g(x-a-|)] (a < x < L) , (12.26) 

sec ^ cos [p (x - L - I)] {L <x < L + a) 

^+(a; + L) = ^+(x) {a<x). (12.27) 

(We may multiply (12.26) by a positive constant factor, and this does not change the 
final result.) We do not write here the expression for x < 0, but only remark that 
ipQ^x) — )■ +00 as X — 7- — oo. On account of the symmetry, i/jq is obtained from i/jq as 

^Q{x)=^+{a-x). (12.28) 

The graphs of V^{x) and V"(x) defined by (11-2) are shown in figure 5. Note that 
V~^{x) — — OO as X — 7- — oo and y~(x) — t- — oo as x — t- +oo. 

Let us calculate go{x, y) and gi{x, y) for < y < x < a. We have 

i'oix) = cosh [s(x — a)] — ^ sinh [s(x — a)], i^oiv) = coshsy + ^ sinhsy. (12.29) 

The Wronskian can be calculated as 

W[iIj^, tpo] = 2^s cosh sa + {l + ^2)s sinh sa. (12.30) 



Low-energy expansion formula for asymptotically periodic potentials 26 

Substituting (12.29) and (12.30) into (11.9) gives go{x,y). To calculate gi, we need to 
evaluate V^"^. It is obvious that V^"^ = Vq~. By definition, 

Q-V+ _ g-Vo" = ^M/P, P= e^p^(^) d^, M= e-^p^(") d^. (12.31) 

J X J X 

Since V^{z) = V^{z) for z > a, we can calculate, by using (12.26), 

/■L+a 111 

P— I r -I-/ MO = - sinhfcg H — sinap, (12.32a) 

J a m{z)? q p 

M = r^\i^t{^)? = ^sech^ ^ + ^ scc2 ^ + i tanh ^ + i tan ^. (12.32&) 

For 0<t/<a;<a, we can also calculate the integrals 

^ v'+(z)j sinh[s(x — a)] sinh[s(7/ — a)] /""^ sinhsx sinhsy 



y si)l{x) sip+{y) ' 5-00 (a;) sV'ol?/)' 

(12.33) 

Substituting the above expressions into (11.86), we obtain 

1 / (V2)sech^ I + (a/2) sec^ f + (1/g) tanh | + jl/p) tan f 
2Y (l/g) sinhfeg + (l/p) sinap 

sinh[s(x — a)] sinh[s(|/ — a)] sinhsa; sinhsy 
5-0+ (x) s^(]"(2/) s'0o'(a^) 



(12.34) 

where i/j^ix), i/joiy) and W^fV'o", V'o'] are given by (12.29) and (12.30). 

The exact Green function for this potential is shown in appendix C (equation (C.2)). 
Comparing with this exact expression, we can check that (12.34) gives the correct first- 
order coefficient of the expansion (see figure 6) . 

Appendix A. Derivation of (7.8) 

We will show that 

/X 
Qix, z\ W: k) gr,(z,u)(x, z: W: k)) dz 
-00 

= lim / Qp{x,z;W;k) gp{z,Up{x,z;W;k))dz, (A.l) 

where gp{x, W) is a function satisfying (4.10). We assume that gp{x, W) is analytic with 
respect to W on the real axis. 
Let a and ^ he defined by 

a{x,z;W;k) /3{x, z;W; -k) 
^{x,z;W;k) a{x, z]W; —k) 

cosh — sinh ^'^^^^^^ \ / a{x,z;k) (3{x,z;—k) 



sinh^^^^ cosh^^^^ / I I3{x,z;k) a{x,z;-k) 



(A.2) 
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Figure 6. The real and imaginary parts of tire exact Gs{x,y;k) for example 2 
(equation (C.2) of appendix C), plotted as functions of real k, with x = 0.4 and 
y = 0.1. All other parameters are the same as in figure 5. The range of k shown in 
these graphs is entirely included in the lowest energy band. (The edge of the band is 
at fc ~ 3.094.) The broken line is the line with slope gi ~ —21.85, which is the value 
given by (12.34). (The value of go obtained from (11.9) is go ~ —3.28.) 



Then, r, Rr and Ri are expressed in terms of a and (3 as 

a[x, z;W;k) 

Substituting (A. 3) into (7.2) and (3.6), we can write Q and u as 

0(x z- W- k) = = ^ = 

' ' ' ^ [a{x,z;W;k) + f3{x,z;W;-k)][a{x,z;W;k) - f3{x,z;W;-k)y 

-( w 7^ T/^ \ , 1 a{x,z]W]k) - (3{x,z;W]-k) 

u{x,z-W-k) = V{z) + log— „. , . ^ J,, — — TT- (A.4) 

a[x, z\ W; k) + p(x, z; W; —k) 

Let a-p and /3p be the quantities obtained from a and /3, respectively, by replacing V 
with Vp. Then, just like (A. 3) and (A.4), we have 

ap[x, z; W; k) 

and 

Qp{x,z-W]k) ^ 



[ap(x, z] W] k) + /3p(x, z] W; -k)\[a^{x, z\ W\ k) - (3p{x, z; W; -k)] ' 

-f i\ T/^ ^ , 1 ap{x,z;W;k) - /3p{x,z;W;-k) 

uJx, z;W; k) = VJz) + log^-- ^ -. A. 6 

' ' ^ ' ^ ap{x,z;W;k)+(3p{x,z;W;-k) ^ ' 

Before carrying out the calculation, let us note that the following equations hold 
for an arbitrary finite number Xq < x: 

lim / Qp[x, z]W]k)gp{z,Cijp{x, z]W]k)) dz 
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fXQ _ nx 

Q^{xQ,z]W]k)g^{z,u)^{xQ,z]W]k))dz+ I g^{z,W)dz, 

(A.7) 



= lim 

fe-)-0 



/X 
Q{x, z; W; k)g-p{z, uj{x, z; W; k)) dz 
-oo 

XO 



/XQ _ pX 
Q{xo,z;W;k)gp{z,u{xo,z;W;k))dz+ / gp{z,W)dz. 
-oo Jxo 

(A.8) 

Equation (A.7) is derived from the two equations 

px px 

lim / Qp{x,z;W;k)gp{z,Up{x,z;W;k))dz^ / gp{z,W)dz (A.9) 

''^^Jxo Jxo 



and 

/xo 
Qp{x, z; W; k)gp{z, Up{x, z; W; k)) dz 
-oo 

/•X'o 
Qp{xo, z; W; k)gp{z, u}p{xo, z; W; k)) dz. (A.IO) 
-co 

The proof of (A.9) is easy. Since the hmit and the integral on the left-hand side are 
interchangeable, we can use (7.7) to obtain (A.9). The outline of the proof of (A.IO) is 
as follows. Prom U{x,XQ]k)U{xo, z]k) — U{x,z]k) and (A. 2), we have 

Q;p(a;, z; W; k) ± /3p{x, z; W; —k) = dp{x, Xq; W; k)[ap{xQ, z; k) ± Pp{xQ, z; —k)] 

+(3p{x,Xo; W; -k)[f3p{xo, z; k) ± ap{xQ,Z] -k)]. (A.ll) 

From (2.1) and (A. 2), we can easily see that 

ap{x, Xo] W] k) ^ cosh ^ ~ ^P^^"-* + Ci{k), Ci{k)^0{k) as A; ^ 0, 

Pp{x, Xo] W; -k) ^ - sinh ^ ~ ^^^"^"^ + C2{k), C2{k) ^ 0{k) as k ^ 0. (A.12) 

We substitute (A. 11) and (A.12) into (A. 6), and calculate the left-hand side of (A.IO) 
by using the method described in section 5 of [2]. Then we find that the contributions 
from Ci(A;) and C2{k) vanish in the limit k ^ 0. It is easy to show that Qp{x, z) and 
u}p{x,z) become Qp{xQ,z) and u}p{xQ,z), respectively, when Ci{k) and C2{k) of (A.12) 
are set to be zero. Hence we obtain (A.IO). Adding both sides of (A.9) and (A.IO) 
yields (A.7). Equation (A.8) can be proved in essentially the same way. Using (A.7) 
and (A.8), equation (A.l) is reduced to 

/xo 
Q{xo, z; W; k) gp{z, lo{xo, z; W; k)) dz 
-oo 

/Xo 
Qp{xo, z; W; k) gp{z, Op{xo, z; W; k)) dz, (A.13) 
-oo 

where xq is arbitrary. The meaning of the equivalence between (A.l) and (A.13) is easy 
to understand. Since limjfc_^oO — linife->oOp hmfc^o*^ — ^^^k^o^p, the difference 
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between the left-hand and right-hand sides of (A.l) in the hmit /c — )■ is determined 
only by the behavior at 2; — > —00. So, in (A.l), we can arbitrarily change the value of 
X as long as it remains finite. Also note that equations (A. 7) and (A. 8), respectively, 
correspond to 

{A;'g^){x,W) = {A;'gp){xo,W)+ r g^{z,W)dz, 

Jxo 

{A-'gp){x, W) = {A-'gp){xo, W)+ I g^{z, W) dz (A.14) 

J Xo 

(see equation (5.34) of [2]). 

Since (A.l) and (A. 13) are equivalent, we will deal with (A. 13) instead of (A.l). 
By taking —xq to be sufficiently large, we can make |Va(-2;)| to be arbitrarily small for 

Z < Xq. 

Let us review how to calculate the right-hand side of (A. 13). (See [2] for details.) 
We split the integral into unit periods as 

Qp{xo, z;W;k) gp{z,Up{xo, z;W;k)) dz = ^A„(/c), (A.15) 
-°° n=0 

rXQ—nL 

An{k) = / Qp {xo, z] W; k) gp{z, ujp{xo, z; W; k)) dz. (A.16) 

J xo — {n+\)L 

Let A and be the eigenvalues of the matrix 



(A.17) 



Q;p(a;, X — L; /c) /3p{x,x — L; —k) 
/3p{x,x — L; k) ap{x,x — L; —k) 

such that |A| > 1 for Imk > 0. Namely, 

A = y - iVl - Y=^[ap{x,x-L;k) + ap{x,x-L;-k)]. (A.18) 
We define 

7 = A-^ (A.19) 

It turns out that A„(A;) depends on n only through 7"'. Moreover, it also turns out that 
An{}i) has the form 

A(A;) = A;7"C(7^A;), (A.20) 

where (7(7", k) is an analytic function of k and 7" on the real axis. We have I7I < 1 for 
Im/c > and I7I = 1 for Im A; = 0. The behavior of 7 for small k is 

7=1 + 2iLoA; + 0{k^) as A; ^ 0. (A.21) 

Since the approach of 7" to 1 is not uniform in n, we cannot change the order of the 

limit and the sum as 

00 00 

hm Vfc7'^C(7",^) = ^lim^7'^C(7",^) = 0. (A.22) 

n=0 n=0 



Low-energy expansion formula for asymptotically periodic potentials 



30 



However, if we treat k and 7" separately, we can let A; — > inside the sum as 
00 00 
^mkJ2rC{r,k) = hmkJ2rC{r,0). (A.23) 

n=0 n=0 

As a result, we have 

00 „i 

lim k y 7"C(7", k)^ / C{x, 0) dx. (A.24) 

fc^o ^ 2Lo Jo 

This equation can be intuitively understood as follows. From (A. 21), we can see that 
y- behaves like exp(2iLoA;n) when k is small. Therefore, 

coo 



poo • nl 

lim k y 7"C(7", k) = lim A; / e^'^°'"'C{e'^'^'''"' , k)diy^— lim / C{x, k) dx. 
fc^o ^ k=o Jo 2Lo fc^o Jo 

(A.25) 



In (A.25), we changed the sum over discrete n to an integral over continuous u, and 
then changed the variable of integration to a; = e'^^^oku^ Here, let us comment about 
the meaning of limfc_^o- The limit A; — ?■ can be taken in any way in the closed upper 
half plane (Imk > 0), but we need to be careful when the limit is taken along the 
real axis (Imk — 0). Recall that we have defined the Green function for ImA; = as 
Gs{x, y; k) = limg^o Gs{x, y;k + ie). So when we write an expression like limfc_^o it 
is to be understood as limfe_>o lim^o (i{k + ie) for real k. Therefore, in (A.25), x — >■ as 
— > 00 inside the limit sign. The limit and the integral in the last expression of (A.25) 
can be interchanged, since \C{x,k)\ is uniformly bounded when |A;| is sufficiently small. 
Hence we have (A.24). 

The left-hand side of (A. 13) is calculated in almost the same way. The only 
difference is that the quantity corresponding to An{k) contains a part which depends 
explicitly on n, in addition to the part which depends on n through 7". If we write 

Q{xo, z- W; k) g^iz, uixo, z; W; k)) dz = J2 Ki^), (A.26) 

then A'^{k) has the form 

A'M = k^Cir, k) + D^ir, k). (A.27) 

As we will see, Dn decreases fike Va{xo — nL) as n — >■ 00. If Va e there exists 

a /c-independent (and 7"-independent) sequence {E^} such that |£)„(7",A;)| < and 
Y^^=o En < 00. Then, we have 

00 00 

hmT,Dn{r,k) = J2Dn{l,0). (A.28) 

n=0 n=0 

We will show that the right-hand side of (A.28) is zero. 

Now let us expficitly calculate the difference between the right-hand and left-hand 
sides of (A. 13). Since |Va| can be made as small as we like by taking large |xo|, we 
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consider only the terms of first order in Va- We can write 

Q{xo, z] W] k)gp{z, u}{xo, z; W; k)) dz 

■oo 

/Xo 
Qp{xo, z; W; k)gp{z, Up{xo, z; W; k)) dz 
■oo 

/xo 
Qa{xo, z; W; k)gp{z, Up{xo, z; W; k)) dz 
-oo 

/xo 
Qp{xo, z; W; k)OA{xo, z; W; k)g'^{z, u}p{xo, z\ W; k)) dz, (A.29) 
■oo 



where 



We define 



QA^g-Qp, UA = u-u:p, g'^{z,W) = ^gp{z,W). (A.30) 



C,p{z) = ap{xo, z; W; k) + /3p(a;o, W\ -k), 

r]p{z) = ap{xo, z; W; k) - /3p{xo, z; W; -k) (A.31) 



and 



C(^) = a{xo, z; W; k) + ^{xq, z; W; -k), 

ri{z) = a{xo, z; W; k) - P{xo, z; W; -k). (A.32) 
Then, as can be seen from (A. 4) and (A. 6), 

Qp{xo, z; W; k) = -, up{xo, z; W; k) = Vp{z) + log (A.33) 

Cp{z)r]p{Z) Cp(^) 

Q{xo, z; W; k) = \ u{x^, z- W- k) = V{z) + log (A.34) 

We can express Qa and ui^ as 

Qa = -7 — (^-\ , = Va - — H , (A. 35) 

CpVp vCp VpJ Cp Vp 

where 

Ca(^) = C(^) - Cp(^), Va{z) = v{z) - r)p{z). (A.36) 
Substituting (A. 35) into (A.29) gives 

/xo 
Q{xo, z; W; k)gp{z,u}{xo, z; W; k)) dz 
-00 

Qp{xo, z] W] k)gp{z, CUp{xo, z; W; k)) dz 

00 

Q{z)lp{z) [^P^^'^^p) +^P^^''^p)] 



rjp{z) = e-"(^) (cosh w A" + sinh w A"") - ike'''^'^ (cosh w^- sinh w A"") p[+]^° 



Low-energy expansion formula for asymptotically periodic potentials 32 

As in (A. 15) and (A. 16), we calculate each integral for xq — nL > z > Xq — {n + 1)L, 
and then take the sum over n. So wc need the expressions of Cpiz)^ Vpi^)^ Ca{z) and 
Va{z) for Xq — nL > z > xq — {n + 1)L. It is shown in [2] that 

Cp{z) = e"(^) (cosh w A" - sinh w A"") - iA;e^(^) (cosh w A" + sinh w A"") pHT'""^ 
+e"C(A)A; + o(A;), 

nL 

+e-''^'^C{X)k + o{k), (A.38) 
where u and are defined by 

In (A.38), and in equations (A.40), we let C(A) stand for a ^^-independent function of A 
which is not necessarily the same everywhere. The expressions for Ca{z) and r]A{z) can 
be obtained, after some calculation, as 

n— 1 

Ca{z) = e"(") 5m(sinhw A"-^"*-^ - cosh w A'^+^'^+i) 

m=0 

n-1 

-ike"^'^ p[-]r""^ Yl ^rnisinh w X^-^^-^ + cosh w X-^+^^+^) 

m=0 

- ^ Va (xq) e"^^^ (sinh w - cosh w A~") 

+ -VA(a;o) e''^^) p[-]^°""^ (sinh ^i; A" + cosh w A"") 

+\ \yA{z) - Va{xq - nL)] e"(^)(cosh w A" - sinh w A"") 
+e"(^)C(A)A; + o(A;), 

n-l 

r7A(^) = e-"(^) Y ^m(sinhw A"-2"^-i + cosh w A'^+^'^+i) 

m=0 

n— 1 

-iA;e-"(^) p[+]r~"-^ XI 5^(sinh w A"-^'"-! - cosh w A-"+2"^+^) 

m=0 

- ^ (a^o) e-"(^^ (sinh w A" + cosh w A"") 

+fv4..)e-«<M+]r-(smh»A"-cosh»A-") 

-J [I^a(-z) - Va{xq - nL)\ e-"(^)(cosh w A" + sinh w A"") 

+e-"(^^C(A)A; + o(A;), (A.40) 



where 



5m = 2 [^a(2:o - mL) - Va{x^ - (m + 1)L)] . (A.41) 
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Let us consider the first term on the right-hand side of (A. 37). We divide the 
integral into unit periods and write 



" -Ca(£) , _ . . ~ r"-"^ -Ca(^) , _ , . 

°° /•XQ—nL 

J2 / a{z,k)b{z,k)d{z,k)dz, (A.42) 

— n J xn — (n+l)L 



n=0 '^xo-{n+l)L 

where we have defined 

, e"(^U"coshw ,/ ,x / _ X ,/ , N -Ca(-2;) 

"^''^"^1(5^' ^(^'^^-^P^^'^^P)' ^(^'^)" e-(-)A--cosh^ - ^^-^'^ 

Prom (A. 38) and (A. 40), we obtain the expressions for a{z,k), h{z,k) and d{z,k) as 

a{z, k) = ao(7") + kai{z, 7") + o(A;), 
6(z,A;) = bo{z) + kbi{z,Y') + o{k), 

d{z, k) = do{z, 7) + kdi{z, 7) + o{k), (A.44) 



(coshw)2(l - co7")(l + Co7")2 

ai(^,7") =iao(7") f c^"^-^ pH^-^ + e"^"^^^ p[+]r"^) +C{r), (A.45) 
V 1 + co7" 1 - co7" / 

^'o(^) ^gp{z,W), 

(A.46) 

n— 1 

7) = 5^ (co7" + 7"-"^) - 2^a(xo) (Co + 7") 

m=0 

[Va{z) - Va{xo - uL)] (1 + C07") , 
d,{z, 7) = -ie^«p[-]^«-"^ ^ (co7" - 7'^-'") 

m=0 

+ '^Va{xo) e^°p[-]^«-'^^ (Co - 7^^) + C(7), (A.47) 

where 

Co = - tanh'iw. (A.48) 

In the above expressions, C(7") and C{'y) stand for the terms independent of z, as in 
equations (A. 38) and (A.40). 

As explained before, we can safely ignore the o{k) terms in (A.44). For any fixed 
finite number /, we can expand 7' in terms of k and neglect the higher-order terms, but 
we cannot do so for 7". We must leave 7"^ as it is. We can see that ao, o-i and 61 depend 
on 7 only through 7". However, this is not the case for do and di. As shown in (A.47), 
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they have the terms including X]m=o ^ml^ and ^^^q -6^7"""^ • We need to know how 
to deal with these terms. The conclusion is that 7*" and 7^™ in these terms can be 
replaced by 1. Let us explain why this is so. Substituting (A. 44) with (A. 47) into the 
last expression of (A. 42), we have the terms involving 7™ and 7""^ as 

00 n— 1 

Y^Sr^ik) J2 Bmicol"^ + 7""™), (A.49) 

n=0 m=0 

where 

Sn{k) = / a{z, k)h{z, k) (1 - i/te-^%[-]^"-"^) dz. (A.50) 

J xo — (n+\)L 

On the right-hand side of (A.50), the terms of order kP vanish since J ao(7")6o(-2^) dz = 
on account of (4.10). So, Sn{k) has the form 

>5„(A;) = A;7"M7"), (A.51) 

where h is an analytic function which can be expanded as h{x) = /io + h^x + h2x'^ + • • •. 
(Note that a{z, k) has a factor 7", as shown in (A. 45). Hence comes the factor 7" in 
front of /i(7") in (A.51).) We can show that 

00 n— 1 00 n— 1 

hm k V 7"/i(7") E = lim ^ V 7"/i(7") E (A-^^) 

n=0 m=0 n=0 m=0 

00 n— 1 00 n— 1 

hm A; J] 7"/i(7") J] 5^7"- = lim k ^Kr) E (^-^S) 

n=0 ?Tt=0 n=0 m=0 

Equation (A. 52) is proved as follows. Since we can express h{'y"') as a power series in 
terms of 7", it is sufficient to consider the case h{'y"') — 7'" with an integer I > 0. We 
need to show that 

00 n—l 00 n— 1 

hm k y 7"('+i) y Bm^"" = lim y 7"('+i) y B^. (AM) 

n=0 m=0 n=0 m=0 

We rewrite the left-hand side of (A. 54) as 

00 n— 1 00 00 

lim A; y 7"('+i) y 5^7'" = lim A; y 5„7'" y 7"('+^). (A.55) 

n=0 m=0 m=0 n=m+l 

Since I7I < 1 for Im/c > 0, we can calculate the sum over n on the right-hand side. (As 
mentioned before, when the limit A; — > of a function a{k) is taken along the real axis, 
it should be understood as limjt_^o o(A;) = limj^^o lime4.o a(A; -|- ie). So we may assume 
I7I < 1 before taking the limit, even when we are considering real k.) We obtain 

lim A; E Bml- E ^"'^'^ = ITo^T^vTT E (^-56) 

m= 

From (A.21) we have 



_ ^ 

m=0 n=m+l m=0 



lim A;^— — = — — . (A.57) 

fe^o l-7'+i 2 / + I L0 ^ ^ 
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Since |7"^('+2)| < 1 for Im/c > 0, and since Y^^=o \Bm\ < oo, we have 



lim V 



m{l+l) 
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(A.58) 



m=0 



m=0 



Therefore, 



n— 1 



2{l + 



n=0 m=0 ^ ' m=0 

It is obvious that the right-hand side of (A. 59) is equal to the right-hand side of (A. 54). 
Hence, we obtain (A. 54). 

It is also possible to derive (A. 52) by using the same argument as in (A. 25). Let 
us define 



B{iJ,) = Bm for m < n < m + 1, 



• /■exp(2iLofei') 



2L 



Then, 



jxi n-1 „oo 

hm k V 7"/i(7") V 5^7"^ = lim A; / e'''^°'"'h{e 

h — ^0 h — ^0 / n 

n=0 m=0 ^ 



Jo 



2iLQki'\ 



m=0 

lim 



h{x) dx. 



(A.60) 
(A.61) 



dFjku) 
du 



Um\F{ku) I B{n)e^'^°'"' dii 



i/=0 



du 



-lim / F(ku)B(u)e^'^'''"' du. 

(A.62) 



It can be seen from (A.61) that 

F(oo) = 0, F(0) 



h{x) dx. 



(A.63) 



2-^0 ^0 

(When k is real, k is replaced by A; -|- ie with positive infinitesimal e.) So the first term 
on the last hne of (A.62) vanishes. In the second term, the hmit and the integral can be 
interchanged since \F{ku)e^^'"°'^'^\ is uniformly bounded and |-6(i^)| du < oo. Hence, 

oo n-1 poo 

lim k V 7"/i(7") y 5^7"* = - / lim F(ku)B(u)e^'^°'"' du 

k-^O ^ — ^ ^ — ^ /n k^O 

n=0 m=0 " 



B{v)du^— / h{x)dxY,Br, 

° m=0 



(A.64) 



This is the same result as (A. 59). Since the last expression of (A.64) is obviously equal 
to the right-hand side of (A. 52), we can see that (A. 52) holds. Equation (A. 53) can be 
proved in the same way. We have 



lim k y 7"('+^) y 5^7"-"* = hm k-^-—. V B, 



/_j_2 ^ • ^ 

7 V- D m{l+l) ^ \ V- ^ 

ml 0/7 9^ r 



n=0 



m=0 



m=0 



2(/ + 2)Lo 



m=0 



(A.65) 
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and 



lim k J2 rnr) E ^^(^) E 

~* n=0 m=0 ° ''O m=0 

oo n— 1 

= lim k E t^/^It'^) E (^-66) 



n=0 m=0 

Thus, we obtain 

oo n— 1 oo n— 1 

lim E ^n(^) E Bmicol"^ + l") = lim E Sn{k) T ^-(co + 7")- (A.67) 

n=0 m=0 n=0 m=0 

Namely, 7"* and 7"™ in (A. 49) can be replaced by 1. 

Now we know that we can replace every 7 by 1, as long as we keep 7" aside. Let 
d{z, k) be the quantity obtained from d{z, k) by this replacement. Then, 

^ rXQ—nL ^ pxo—nL 

lim / a{z,k)b{z,k) d{z,k) dz = lim / a{z,k)b{z,k) d{z,k) dz. 

(A.68) 

With 7=*="* — >■ 1, we have 

n—l n—1 ^ 

E Bml"" ^ E = 2^^''^''°^ ~ ^""^^^ ~ 

m=0 m=0 

n—l Ji— 1 ^ 

E Bmr-"" ^ E = 2 ~ ~ ^"^-^^^ 

m=0 m=0 

and so equations (A. 47) become 

do{zn) ^ 4(^,7") + ^0(7"), d,{z,^) ^ Ji(z,7") + C(7"), (A.70) 

where 

dA{z, 7") = -^Va(^)(1 + C07"), (io(7") = ^Va(xo - nL)(l - co)(l - 7"), 

d,{z, 7") = '-Va{xo - nL) e^Opf-]^"-"^ (cq - 7") , (A.71) 

and C(7") is an ^-independent term. With (A.71), we can write d{z, k) as 

d{z, k) = dAiz, 7") + + kdi{z, 7") + A;C(7"). (A.72) 

We substitute (A.44) and (A.72) (with (A.45), (A.46) and (A.71)) into the right- 
hand side of (A.68), and calculate the integral over z. Since gp satisfies (4.10), we 
have 

fXQ—nL fXo—nL 



/ ao(7")&o(^)fio(7") = a,{^^)d,{^^) / g^iz, W) dz = 0, 

J XQ — (n+l)L J xo — {n+l)L 

/ aombo{z)C{r)dz^ao{r)C{r) / gp{z,W)dz^O. (A.73) 
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Therefore, 

/ a{z,k)biz,k)d{z,k)dz = Fr,{-f'') + kGni^'',k), (A.74) 

J xo—('n+l)L 



'xo—{n+l)L 

where we have defined 

pxo—nL 

F,(7") = aoir) / bo{z) dA{z, 7") dz, (A.75) 

"XQ—nL fXo—nL 



rXQ—riL, pxo—nij 

Gr,{r,k) = do{r) / ar{z,r)bo{z)dz + ao{r) / d,{z,r)bo{^)dz 

J XQ — {n+l)L J X{) — {n+l)L 

pxa—nL 

+ao(7") / dA{z,r)bi{z,r)dz + 0{k). (A.76) 

Note that F„(7"') + A;G„(7", k) is a part of Dn{'y'^, k) of equation (A. 27). Since c^a, c^o 
and include either V/^{z) or VA(a;o — nL), both and Gn decrease hke V/^{xq — riL) 
as n — )■ oo. (The 0{k) term in (A.76) also decreases like VA{xo — nL).) So, if Va G Fq~\ 
then, as in (A. 28), 

oo oo oo 

^^Y}F^m + kGn{l\k)] = J]lim[F„(7'^) + A;G„(7^A;)] = 5]F„(1), (A.77) 

n=0 n=0 n=0 

where we have used |Gh(1,0)| < oo. Substituting 

ao(l) = d^iz, 1) = -^(1 + c,)Va{z) (A.78) 

i + Co z 

and hQ{z) = gp{z, W) in (A.75), we obtain 

Fnil) = -t: / ^a(^)^p(^, W) dz, (A.79) 

^ Ja:o-(n+l)L 

and hence 

hmV / aiz,k)b{z,k)d{z,k)dz = y2F^{l)^-- / V^a(-2)5p(^, W^) d^. 

(A.80) 

Thus we have obtained, for the first term on the right-hand side of (A. 37), 
Since gp{x, W) also satisfies (4.10), equation (A. 81) holds with gp replaced by g'^: 
In the same way, we can derive 
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It is easy to see that 

lim / —— - — —gAz,u:^)az= / lim — — - — —gAz,u:y)az 

yA(^)^;(^,W^)dz. (A.85) 



Taking the hmit /c — )■ of (A. 37), and substituting (A.81)-(A.85), we obtain (A. 13). 
This conclusion does not change when the terms of higher order in Va are taken into 
account. 

Appendix B. Proof of (7.15) 

Since U {z, x; k) is the inverse of U (x, Z] k), ioi z < x we have a{z, x] ±/c) = a{x, z\ ^k) 
and (3{z,x;±k) = —(5{x, z;±k). Using this, from (2.1) we obtain 

d _ 

—a{x, z; W; k) — ika{x, z; W; k) — f{z)P{x, z; W; —k), 

^^(x, z; W; -k) = -ikp{x, z; W; -k) - f{z)a{x, z; W; k), (B.l) 
where a and P are defined by (A. 2) of appendix A. Hence, for Im/c > 0, 

^{\a{x,z;W;k)\^ -\f3{x,z;W;-k)\^) 

= -2(Im A;) {\a{x, z; W; k)\'^ + z; W; -k)f) < 0. (B.2) 
Since \a{x,x;W;k)\'^ - \p{x,x;W; -k)\'^ = 1, from (B.2) it follows that 

\a(x, z; W; k)\'^ - \p(x, z; W; -k)\'^ > 1 (B.3) 
for z < X. Using (A. 3) of appendix A, we have 

\f{x, z; W- k)\^ + \Ri{x, z- W- k)\^ < 1, (B.4) 



and this gives 



\Q{x,z;W;k)\^ 



f\x,z;W;k) 



\r{x, z; W; k)\ 



2 



< — '^1' ' ' ^' < 1. (B.5) 

- l-\Riix,z;W;k)\^ - ^ ^ 



1- R^{x,z]W]k) 

Appendix C. Exact Green function for Vs given by (12.23) 

Let us define 

Pk = + Eo, Qk = \/C - k"^ - Eo, Sk = \/ h - k"^ - Eq, 

a = e^''^ I cospfeO + — — — sinp^a ) , a' = e"^*^** ( cospfcO - — — — sin pka 
\ ^PkQk J V ^PkQk 

Ci = - e '^''"sm Pka, C2 = ^ ^ -. C.l 

The exact Gs{x, y; k) for (12.23) if) < y < x < a) can be obtained as 

Gs{x,y^k)^- (C.2) 
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where 



(Ci + C2) cosh[sfc(a; - a)] + {qk/sk){Ci - C2) smh[sjk(x - a)], 
(Ci + C2) cosh SfeX - {qk/sk)iCi - C2) sinh SfeX, 



(C.3) 



= [(Ci + C^f + (gfcM)'(Ci - Sksinska - 2qu{Cl - Cl) cos s^a. 

(C.4) 
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